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Abstract 

We study the spin fluctuation of the oxygen deficient planes of YBa2Cu306+a;- 
The Cu-0 chains that constitute these planes are described by a model that 
includes antiferromagnetic interactions between spins and Kondo-like scatter- 
ing of oxygen holes. The spectrum of magnetic excitations shows the presence 
of incommensurate dynamic fluctuations along the direction of the chains. 
The presence of itinerant holes is responsible for the existence of important 
differences between the spin dynamics of this system and that of a quasi- 
one-dimensional localized antiferromagnet. We comment on the possibility of 

experimental observation of these fluctuations. 

PACS numbers: 75.10.Jm, 75.40. Gb, 74.72.Bk, 61.12.Ex 
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I. INTRODUCTION 



A considerable amount of experimental and theoretical research has recently been con- 
secrated to the study of the structural and electronic properties of the oxygen deficient (or 
CuO^) planes of superconductors of the YBa2Cu306+x-type. It has been known for some 
time that these planes play an important role in the superconductivity of the YBCO com- 
pounds by regulating the hole concentration in the Cu02 planes There is now evidence 
that they also have an interesting physics of their own The basic components of the 

CnOx planes are Cu„+iO„ chain fragments in which Cu and O atoms form a linear alter- 
nating chain. Under stoichiometric conditions, x = 0.5 and x = 1, all chains are ( ideally) 
infinitely long. Away from stoichiometry the length of the chains follows a certain distribu- 
tion determined by x and the temperature but the physical properties of the chain planes 
are dominated by the longest chain fragments Band structure calculations for 

YBa2Cu307 show that the dispersion of the chain-derived bands which cross the Fermi level 
in the directions perpendicular to the chains is much weaker than that in the direction of the 
chains. This suggests that, to a first approximation, the latter may be viewed as decoupled 
from each other and from the Cu02 planes. If this is correct, important 1 — d fluctuation 
effects are expected in the charge and spin excitation spectra of the CuOa, planes. A recent 
inelastic neutron scattering study of a YBa2Cu306.93 crystal shows that such effects are 
indeed visible in the phonon spectrum of the chain planes. Their signature is the presence 
of purely inelastic satellites around the structural Bragg peaks at positions that correspond 
closely to the value of 2kF for the one- dimensional chains predicted by the band structure 
calculations 0. 

In this work we study the spin fluctuation spectrum of long Cu„+iO„ chains in the 
framework of a one-dimensional model of strongly interacting fermions. Most of the paper 
deals with the T = case but we briefly discuss the extension of our calculations to finite 
temperatures in the concluding Section. Our main results are the following. The spin-spin 
correlation function of the chains is given by the convolution of two functions. One is the 
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dynamic structure factor of a one-dimensional Heisenberg antiferromagnet on a fictitious 
lattice that contains a number of sites equal to the total number of holes in the chain. This 
factor describes the part of the dynamic that results from the exchange interactions between 
the spins in the chain. The second factor, related to the density- density correlation function 
of the itinerant oxygen holes, takes into account the perturbation of the spin configurations 
that arise from the motion of the latter. The magnetic intensity is modulated along the 
direction of the chains and is peaked at an incommensurate position, Qmax = vr(l + c), 
where c is the nominal oxygen hole concentration in the chains (see below). The lineshape 
of the energy-integrated peak, very asymmetric, is quite different from that obtained for a 
one-dimensional Heisenberg antiferromagnet. The intensity at Qmax is distributed over an 
energy range given by the Fermi energy of the holes. However, enough intensity remains at 
low energies that observation of these excitations by inelastic neutron scattering might be 
feasible. 

The organization of the rest of the paper is as follows. The model used to describe the Cu- 
O chains is presented in Section ^ The calculation of the ground state spin-spin correlation 
function is described in |TT|. In Section |V| we discuss our results and their connection with 
experiment. 

II. THE MODEL. 

The simplest microscopic model for the Cu-0 chains is the 1 — d Emery Hamiltonian 
defined by: 

^ = E E {pI^^ + dUPra) + E(e 

pUr + UpHr^rirl) + E (^-i^P + Udrip^Upl) , (1) 
<r,p> a r p 

Here, r and p denote the O and Cu sites, respectively, d„ and p„ are annihilation operators 
for holes of spin a on Cu and O sites, the occupation number of a hole with spin a is and 
n = J2a''^cr- Hopping is only allowed between nearest neighbors. The values of the charge 
transfer gap and of the on-site Coulomb repulsion appropriate for the description of the 
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CuO chains are (e^ — e^) = 2t, and Ud = 2Up = 8t, where t ~ leV is the Cu-0 transfer 
integral. In the following we shall measure energies from the position of the Cu level and 
set ed = 0. 

An additional parameter, the charge of the fragment, has to be fixed in order to com- 
pletely define the problem. Taking the Cu^^ and configurations as reference states, one 
may easily see that a neutral fragment with Cu atoms contains 2N — 1 holes. However, 
because of the transfer of charge between the CuO^^ and the Cu02 planes the number of 
holes in the chains, A^(l + c) with <c < 1, deviates from this value |jlO|. In a naive picture 
in which all the Cu ions in the chain are in the Cu^^ state, c is simply the concentration 
of oxygen holes. It is experimentally known [|ll|] that the hole concentration in long chain 
fragments is u^^LS. Therefore, for the systems discussed here, c~0.3. It is however inter- 
esting to treat c as a free parameter and study the dependence on concentration of different 
physical quantities. 

The Emery Hamiltonian of Eq. |I] gives a good a description of the static electronic 
properties of the chain fragments |TD|,|T^ but it can only be properly treated with numerical 
methods such as Quantum Monte Carlo or exact diagonalization that are not well suited for 
the calculation of the dynamic susceptibility. There is a limiting case, Ua/t — > oo, ep/t and 
Up/t finite but large, in which a more tractable model [|l^] can be derived from Eq. [^. In 



this limit the Cu occupation Up is unity and charge fluctuations on the Cu sites only appear 
in virtual processes generated in perturbation theory by the kinetic energy term. In the 2nd 
order in t we generate a process in which an O hole hops from a site r to an adjacent O site. 



r±2 leaving the spin sequence in the lattice unchanged [Q. This process, whose amplitude 
is T = t'^/ep, is described by a tight-binding model for spinless free fermions defined on the 
oxygen sublattice: 

N 

Hi = -r E(c]+i+i/2S+i/2 + ^-c-)- (2) 
i=i 

In the same order in t there is an exchange interaction between a spin on an oxygen site 
and the spins on neighboring Cu sites with coupling J = t^/ (cp + Up). Finally, if there is 



no O hole on a site r the Cu spins on sites r ± 1 still interact antiferromagnetically with a 
coupling constant that is of the fourth-order in t, J' = At^ /{e^{2ep + Up)). The spin part of 
the interaction is thus described by a 1 — d Heisenberg Hamiltonian given by : 

H2 = ^ Jj^i+lCTj.CTj+i, (3) 

i=l 

where Nh = Nc is the number of oxygen holes in the system, + Nh is the total number of 
spins in the system and cTj are spin 1/2 operators. The nearest-neighbor coupling Jj j+i = J 
or J' depending on whether + 1) is a Cu-0 or a Cu-Cu pair, respectively. 

Since the couplings JiA+i depend upon the instantaneous hole distribution, the spin and 
charge degrees of freedom of the system are coupled and the total Hamiltonian is not exactly 
solvable in general. However, in the special case J' / J = 1 an exact solution can be found. 
This is because if the exchange coupling is the same for all pairs of spins regardless of the 
charge configuration, the Hamiltonian becomes the sum of two independent terms each of 
which is exactly solvable. In the case of the Cu-0 chains the ratio of couplings J' / J ~ 0.75 is 



not unity but, as shown elsewhere |T0|, their statzc magnetic correlations are quite accurately 
described by an effective homogeneous model. It seems therefore natural to try and use it 
to compute the time- dependent correlation functions as well. This is the approach that we 
follow in this paper. 

The appropriate parametrization of the effective model may be determined variationally. 
If we construct a trial ground-state wave function as the product of separate wavefunctions 
for the spins and for the holes and then determine the individual components by minimizing 
the total energy of the system, we obtain the effective model 

N _ N+Nh-1 

Heff = E(sVl+l/2Ci+l/2 + ^^-C-) + Jeffic) ^ CTi.CT.+i. (4) 

j=l i=l 

where Jeff{c) = J'{1 — c) + c J and we have omitted an energy shift that can be absorbed in 
a redefinition of the chemical potential of the holes. We find Jg/j/r ^ 0.3 for c = 0.3. 
In spite of the simple form of H^ff the spin-spin correlation function, 

/oo 
rfte^™(SQ(t).S_Q(0)), (5) 
-oo 



is non-trivial because the relationship between Sj and cTi is not simple |jT5|]. The difference 
between the coordinates of a spin measured in the real and in the fictitious lattices is 
determined by the number of holes located to its left. In terms of cr^,, the Fourier transform 
of the spins on the fictitious lattice, we have: 



i-i 



S, = (iV + iV,)-V2 ^ exp ik{l + ^ c]^i/2S-+v2) 

k 3=1 



(6) 



for Cu sites, and 



i-i 



Si+l/2 = (A^ + Nh) exp ik{l + 1 + ^5+1/2^+1/2^ ^/+l/2^«+l/2, 



(7) 



k 3=1 

for O sites that only carry a spin when they are occupied by a hole. 

Substituting Eqs. || and |^ in Eq. |^ and using the fact that the ground-state wavefunction of 
the decoupled Hamiltonian factorizes, S{Q, Vt) may be written as a convolution : 



SiQ,n)= — \A{Q,k)\'f{k,u;)a{Q-k,n-u;) 

-TT ZTT J-oo /TT 



(8) 



where k = k{l + c). 

The first factor in Eq. H, f{k,u!), is the dynamic spin correlation function of the fictitious 
1 — d Heisenberg model of Eq. S, 



f{k,u) 



dte'^'{cTkit).(T^km 



H2 



(9) 



The second factor, C^, describes the effects of the hole dynamics and is given by the space- 
time Fourier transform of : 



Ck{lJ';t,t') = ^exp 



ikJ2<l>{qJ)p{q,t) 



exp 



-zA;^0(g,r)p(g,t') 



(10) 



Hi 



where the density operator is given by : 



(11) 



and we have defined (f){l, g) = (1 — e*''')/ (1 — e*') . The g = component of the density is not 
present in Eq. [TD|. Searate treatment of this component leads to the appearance of k instead 
of k in the argument of Ck in Eq. |. Finally, A{Q, k) is a coherence factor given by 
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This factor takes into account the phase shift between Cu and O spins and it is derived in 
Appendix 0. 

III. EVALUATION OF THE CORRELATION FUNCTION. 

The evaluation of the dynamic structure factor in Eq. ^ requires the knowledge of the 
individual correlation functions f{k,uj) and Ck{q,uj). The dynamic correlation function of 



the one-dimensional Heisenberg antiferromagnet has been discussed in reference |]16|. This 



function cannot be obtained in closed form but a very accurate approximate expression for 



it has been derived [T^ using a combination of exact analytic and numerical results : 

f{k,uj)=3ae{uj)e{7rJ,ffSm{\k\/2)-uj) (u^ - {^J,ff sin kf^ (13) 

where a = 1.4579 is a numerical constant. There are no sharp spinwave excitations in the 
system but the spectral function is strongly peaked at the lower edge of the spectrum where 
the spectral weight diverges as (to — ujs{k))^^^'^ with uJs{k) = 7rJeff/2smk. 

It will be seen below that most of the spectral weight in S{Q, Q) is concentrated in the 
neighborhood of Qmax = 7r(l + c) (mod 2tt). Since f{k,uj) is strongly peaked, the structure 
factor in the neighborhood of Qmax is determined by the long distance part of the fermionic 
correlator Ck in Eq. ||. This part of the response is determined by the low energy excitations 
of the electron gas in which the momenta of excited particles and holes are all near ^kp, 
the Fermi momentum of the holes. 

Eq. may be evaluated using a well known procedure Jl^ in which the fermion density 



is split into contributions from right and left moving fermions with constant group-velocity. 
Since the fermionic Hamiltonian is particle- hole symmetrical it is sufficient to compute Ck 
for c< 1/2 and make the transformation c — (1 — c) to obtain the corresponding expression 
for c > 1/2. We write 
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p>0 p<0 

and approximate the group velocity of right and left moving fermions by ±vf respectively, 
where the Fermi velocity is Vp = 2rsin/ci? and kp = tcc {kp = 7r(l — c) for c > 1/2) . 
With this approximation the time dependence of the density operators becomes very simple, 
P±(q') t) = e^*'^''*p-|- (g) and ujg = vp \q\. 

The long wave-length density fluctuations of the 1 — d electron gas may be represented 
by boson operators defined by the relations |TP[ : 

P+iq) = viq) {o{q)b+ + ei-q)b_;) ^^^^ 

p-{q)=v{q) (9{q)h_g + e{-q)h+), 



where v{q) = J\q\/ {2Ntx). With the definition 



(16) 



we may rewrite Eq. ^ in the form : 



Ck{l,l';t,t') = {expikA{l,t) exp-ikA{l',t'))^^ 
The evaluation of this expression in Appendix yields : 



Ck{l,l';t,t') = exp 



An' 



where 



(17) 



(18) 



\q\<:kp 



-ILOnt I 



1 - COS q{l - I')) + (1 - e"*"^'*) (2 - cos ql - cos qV 



■ (19) 



2(1 — cosg) 

The integration over momenta in Eq. |1^ is cutoff at q kp, the wave vector below which 
the density fluctuations of the electron gas can be treated as well defined excitations. In the 
thermodynamic limit (/,/' — oo, |/ — l'\ finite ) J-'{l,l';t) depends on / and /' only through 
the difference (/ — /') . Evaluating the integral over momentum in Eq. |TU| in the logarithmic 
approximation we find : 



Ck{l; t) = (1 + ikp{l + vpt)){l - ikp{l - vpt)){l + iEptf 



(20) 



where Ep = vpkp is the Fermi energy, and ak = (fc/27r)^ < 1/4. 

The space-time Fourier transform of Cfc(Z; t) can be found analytically. For q in the first 
Brillouin zone : 

A 2n(. .\ (ou \1— 4cn; —u}/Ep 

Ep T{ak)T{3ak) [uj^/vj,-q^Y ■iak\q\2.ak 
where F is the Gauss hypergeometric function. Ck{q,uj) diverges at the lower end of the 
spectrum, at u; = 1^1^^, the dispersion relation for density waves in the electron gas. At high 
energies the spectrum is cut off exponentially at ^ Ep. 



The instantaneous correlation functions can be obtained by integrating Eqs. [T^ and pi 



over frequency. For the Heisenberg model the result is [|16 

1 + sin|A;/2| 



m = ^ log 

Zn 



3a 
Zn 



n — \k\ 



(22) 



cos(A;/2) 

The approximate expression quoted on the right hand side of Eq. ^ derived near the zone 
boundary at k = ±TT where f{k) diverges logarithmically is in fact remarkably accurate over 
most of the Brillouin zone and will be used in our calculations below. The instantaneous 
value of the fermionic correlator is given by : 



Ck{q) 



2^-"* \/2tt 



kp 



K,/2-a,{q/kp). (23) 



T{ak) kp 

Ck{q) is strongly peaked at g = where it diverges as Ck{q) ~ ^^(1^20^)^ Away from the 
center of the Brillouin zone Ck{q) decays exponentially, Cfc(g) ~ exp{—\q\/kp). 

The rapid decrease of f{k,uj) away from the zone boundary restricts the interval of the 
/c— integration in Eq. ^ to the vicinity of A; = ±7r. This allows us to linearize the bottom of 
the spectrum of the spin excitations writing uJs{k) ~ 7rJe//(7r — |^|)/2 and to replace the 
exponent in Eq. ^ by a = 1/4, its value at the zone boundary. With this approximation 
the fermionic correlator becomes : 

^ / N ^/ .4:'Kexp(—uj/Ep) 1 ^/ s 

Carefully analyzing the contribution to the integral over k in Eq. ^of the different regions 
where the spectral weights of / and C are simultaneously large, and using the fact that the 
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structure factor varies smoothly compared to the two other factors in the integral, we may 
write the dynamic structure factor in the approximate form : 



s{Q, n) = |A+|2$(g - 7r(i + c), n) + \a_\^^-q + n{i - c), n), 



(25) 



where we have defined 

duj e"""/-^^ f°° dq 



$(/t,fi) = 3a 



VsEpil + c) 



ixJ2\q\ L 



Vt — uj 



v.. 



n-l/2 



\q — K\ 



— ] -q- 

Vp, 



.-1/4 



• (26) 



Here, Vg = 7rJe///2(l + c) is the renormalized spinwave velocity and the amplitudes A± are 
given by 



TTC 

A± = l± sin y. 



(27) 



Eq. ^ contains the main result of this work. The structure factor of the chains appears 
as the sum of two contributions. We will be seen below that $(K,fi) peaks at k ~ and 
is exponentially small for k > 0. Therefore, the first term is important for Q ^ 7r(l + c) 
and the second one for Q ^ 7r(l — c). For vanishing hole concentration A± = 1, and the 
two $ functions reduce to the Q < n and Q > tt parts of f{Q, fi), respectively. For a finite 
hole concentration the corresponding peaks shift and broaden and the result is a distorted 
lineshape. 

In order to make further progress we must analize in detail the properties of ^{k,Q). 
We begin by examining the properties of the instantaneous function obtained by integrating 
Eq. ^ over frequency. The result can be expressed in terms of the convolution of the two 
static functions defined in Eqs. |22| and B3, 



$(k) 



3a 



dq 



C(g)log 



4(1 + c) 



q — K 



(28) 



27r(l + c) A 27r 

where C{q) = C^(g) (cf. Eq. ^) and we have used the asymptotic form of f{k) given in Eq. 
22. 

For c ^ and c 1 where kp vanishes, C{q) 27T6{q) and $(k) 3a/(27r(l + 
c))^(— k) log |4(1 + c)/k|. For intermediate values of the hole concentration $(k) may be 
obtained analytically in two limits. At large we find : 
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3a 



27r(l + c] 
3akp 



log 



4(1 + c) 



K 



21/4 



^ 167r(l + c)2 r(l/4) Kkp 



-3/4 



-K/kp 



7 - 1 + log 



4(1 + c) 



, K <^ —kp 

,K:§>kF 



(29) 



where 7 is the Euler constant. We thus recover the behavior of the pure system for large 

negative values of n. The exponential tail that we find for large and positive k results from 

the smearing of the unperturbed correlation function by the moving holes. 

We can also evaluate analytically Eq. ^for small n/kp- The result is : 

6a 

I - 7^ 

$(k) 



(1 + 1/2 log|(l + c)^|) 



(27r3)V2(i + c) " ' ' '^^"^ ' '"^^ 
•^(0) - TTT^w^^^ V^^T^ (7r/2-l-l/2 log |(1 + c)/«:|) , < 



(27r3)l/2(l + c) 

where $(0), the value at k = 0, is 

$(0) 



(30) 



3a 



\ 4(1 + c) 
log ^7 + Co 



(31) 



47r(l + c) \ kp 

and Co = 7 + 7r/4 + 3/21og2 ^ 2.4. For a finite concentration of holes $(k) has a cusp at 
K = 0, a less singular behavior than at c = where it diverges logarithmically at k = 0. 

The two asymptotic forms, Eqs. |2S| and match smoothly at k, ^ kp as shown in Fig. 
ID where we plot the k— dependence of $(k) for different hole concentrations. 

We now turn to a discussion of the frequency dependence of the spectral function. The 
excitation spectrum of the pure model goes down to zero-energy only at k = 0. The same is 
true for a finite hole concentration. Putting k = in Eq. ^ and performing the integral we 
find : 

1 



$(0,^) 



1 



2Ef (1 + c) Jo 

r(i/4) 



—— ^(x,fi) F (1/4, 1/2,1, 



X + 



dx 



V^(x,(]) f(i/4,1/4,3/4,/i2(x)), (32) 



2kpy/27rVsVp r(3/4)(l + c) J{i+vs/vf)-^ ^^{1 - x) 
where we have defined the functions ^p{x,^l) = 6(x — 1 + nJeff /^) e~^^^^f, and h{x) = 
{vp/vs){l — x)/x. Fig. shows a plot of $(0,fi) as given by Eq. |3^ for several hole con- 
centrations and for r/J^ff = 3, a value consistent with our parametrization of the Emery 
model. 
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In the pure Heisenberg model the spectral weight, $(0,ri) ^ diverges at zero 

frequency. This divergence is cutoff for finite c by the motion of the holes and the intensity 
at i7 = is now finite. The weight missing at low energies is transferred to the high energy 
end of the spectrum where an exponential tail, $(0, Q) ~exp {-n/EF){nJeff/n)'^/\ extends 
beyond the energy Qmax = ^rJe// where the spectrum of the pure Heisenberg model has a 
sharp cutoff. Analysis of Eq. |26| for finite k, is more complicated. In this case the square 
roots present .in the integrand of Eq. ^ impose the inequalities VfIqI < uj < Q — Vs\q — k,\ 
on the domain of integration. This condition can only be satisfied if i7 > K,mm{vF,Vs), i.e., 
there is a threshold in ^{n,Q). Let us first consider the case vp > Vs- The bottom of the 
spinwave continuum starts aX Q = Vs\k,\- Near the onset, for Af2 = fi — Qmin ^ only the 
neighborhood of g = contributes to the integral in Eq. ^ Making the change of variables 
q = {uj/vf) sin 6' the latter may be rewritten as : 

Jo E/(l + c)27r ,/ tan^ V 



X (^Afi-cu(l-^sine)^ . (33) 

The first factor in the integral over 6 is dominated by 2kVs ^ Q ^ AQ which allows us to 
approximate 

^U,Q) ^ — {AQ-u(l smf 

Jo Ef{l + c)2n^/KVsJ^n/2 itanei^ '"f 



VAn W2 de 

' (34) 



Ef{l + c)tt./M; J-n/2 J\tan9\(l - ^ sin 



Vp 



The last integral can be done and we find the result 

^ 2 r^)"''^-";/"^)-"'. (35) 

^ ' ^ \2kVsJ Ef{l + c) ^ ' 

In the case Vs > vf, the bottom of the excitation spectrum is at Q = vf\k\, i.e., below 
the spinwave continuum. The integral over q in Eq. ^ is now dominated by g ~ k. The 
integration may be performed in a way similar to the previous one and the result is : 



^ ' ' 3 \2kvfJ Ef{l + c) 



-K/kf 



(36) 
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The motion of the holes modifies substantially the low energy part of the excitation spectrum 



of the chains. For finite c, vanishes at threshold (cf. Eqs. ^ and ^) in contrast 

with the situation at c = where the spectral weight has an inverse square root divergence 
at r2 = Vs\k\. For c~0 and c~ 1 where vp < Vs, there are hole-like excitations below the 
spinwave continuum. However, their total weight is small because they are confined to the 
region |g— 7r(l+c)| ^ A^^by the exponential factor in Eq. Elsewhere in momentum space 
the spectrum is essentially that of the Heisenberg antiferromagnetic chain. At intermediate 
hole concentrations Vp > Vg and the excitation spectrum starts at the spinwave continuum 
but there exist excitations with energies up to Ep that drain spectral weight from the region 
uj ^ Jeff- The crossover between these two regimes takes place when Vg = vp, which for our 
parametrization of the model yields a critical concentration Ccrit ~ ^ sin~^(7rJ'/4r) ~ 0.06. 

IV. DISCUSSION 

The natural tool for experimental investigation of the magnetic fiuctuations described 
above is inelastic scattering of spin polarized neutrons. In this Section we will use our 
results to compute the differential cross-section for magnetic neutron scattering. The latter 
is related to the dynamic structure factor of Eq. |] by the expression 

where tq = e^/mc^, and 7 = 1.57 is the neutron's gyromagnetic ratio. The neutron's initial 
and final momenta are ki and kf, respectively, q is the momentum transfer in the direction 
parallel to the chain, u is the neutron 's energy loss and dQg is the element of solid angle in 
the direction of the outgoing neutron. 

The intensity scattered by a quasi one-dimensional object is uniformly distributed in a 
plane perpendicular to it. Therefore, the intensity at any particular point in g— space is 
very small and it is generally impossible to resolve the structure factor both in q and u. 
Under typical experimental conditions what is measured is the cross-section integrated in 
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energy from up to the incident neutron energy, ujn- If the experiment is performed under 
quasielastic conditions (i.e., if most of the spectral weight is at energies smaller than a;jv) 
we have 

^ = \{r,^fS{q) = 194 mbarn 5(g), (38) 

where is has been assumed that j^^ du S{q, u) ~ /q°° du S{q, uj) = S{q). 
Using the results of the preceding Section we may write 

S{q) = W+^{q - 7r(l + c)) + W^^-q + 7r(l - c)). (39) 

where W± = A^, and $(k) was defined in Eq. |28|. Figure ^ shows the evolution of S{q) 
with doping. In the limiting cases c = and c = 1 the system is equivalent to a Heisenberg 
antiferromagnet. In the undoped case, = = 1 and the two terms in Eq. ^ contribute 
to the peak at g = vr, one for q<TT, the other for q>TT (cf. discussion below Eq. |2^). In the 
fully doped case =4, =0 and only the first term remains giving a peak at g = 27r. The 
fourfold increase in the amplitude that occurs when going from c = to c = 1 is halved 
by the factor (1 + c)~^ that appears as a prefactor in the expression for $(/t). Therefore, 
the area Jq^ dq/{2-K) S{q), proportional to the number of spins per unit cell, just doubles 
as it must. As the hole concentration increases from c = to c = 0.5 the peak shifts and 
the lineshape becomes increasingly different from to that at of the Heisenberg chain because 
the interval of momenta in which the presence of holes is relevant, Aq ^ kp, gets bigger. 
Upon further increase of c the process reverts because kp decreases back to zero (cf. the 
paragraph above Eq. |1^). 

The hypothesis of quasielastic scattering implicit in the precceding discussion is unlikely 
to be fulfilled under realistic experimental conditions. For c ~ 0.3, the excitation spectrum 
extends up to energies of the order of Ep ~ leV, whereas u^- ^ 50 meV in a typical 
spin-polarized neutron scattering experiment. An important reduction of intensity is thus 
expected to result from the fact that the excitation spectrum is only partially integrated. 
The peak intensity can be estimated by integrating Eq. ^ at g = 7r(l + c) up to at, using Eq. 
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32| . Wefm(].da/dVLs\max^ 60 mbarn for cjjv:^ 50 mev whereas, under the same conditions, Eq. 
^ would predict da / dVLs[max^ SOOmbarn. However, despite the large decrease of intensity, 
the predicted signal might still be observable in a high flux reactor. 

We conclude by examining the expected finite-temperature corrections to the results 
obtained in this paper. The effects of temperature are different in the fermionic and spin 
subsystems. We show in Appendix |C| that the power-law behavior of the correlation function 
Ck{l',t = 0) that we found in Eq. gives way to exponential decay at finite T : 

Cfc(/;t = 0) = expf- — j , (40) 

where the typical scale of the correlations is ^ = vp/T. The appearance of a finite decay 
length rounds-off the singularities of C{q, u) at a scale q ~ co'~T. At room temperature 
T ~ 300 K, and for f ^ ~ 1 eV this should lead to a smearing of the peaks in Fig. ^ on a 
scale ~ 0.03r.l.u., which is of the order of magnitude of the experimental resolution. 

The changes in the spectral properties of the Heisenberg chain at finite T have been 
discussed by Miiller et al. |T6[ who have shown that spectral weight appears below the 
bottom of the spin-wave continuum with increasing the temperature. At ksT/J = 0.5, 
the total intensity in the region u < Ugik) , is comparable to that aX u > Ugik) but for 
ksT / J ~ 0.1 these effects are negligible and the correlation function is essentially identical 
to that at T = 0. For CuO chains where J^jf — 0.15 eV and if the experiment is performed 
just above the superconducting transition temperature at T 100 K, the ratio is T/ J ~ 0.1. 
Under these conditions the predictions of our T = calculations should remain valid. 
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APPENDIX A: THE COHERENCE FACTOR 



In this Appendix we derive the form of the coherence factor A{Q, k) defined in Eq. 12. 
The Fourier components of the spin density are given by |1^ 



Sq = E e-^^' (S. + e-«/2S,+i/2) = Sq|^„ + SqIo , (Al) 
I 

where the right hand side of the equation defines the Cu and O contributions to Sq. Using 
the expansions ^ and ^ above and the identity h = {expikn — 1) / {expik — 1) vahd for an 
y number operator n and any value of k, the oxygen contribution may be cast in the form : 

SqIo = iN + iV.)-V^ ^ ^ exp ^ [{k - Q){1 + 1/2) + k^{l)] """^^^^^^^J^ (A2) 

where (^{l) = Z]j=o%+i/2- Noticing that Lp{l) + = ^p{l + 1) this is equivalent to 

S« lo = ( + A'.)-"^ E E exp - g) + y^) ^''^^^^/^''^^''' ■ (A3) 

Changing indices of summation from / to /' = / + 1 in the first summation in the equation 
above and recalling that periodic boundary conditions impose that f{N) — (f{0) = Nh and 
that kN{l + c) = 27rmi and QN = 27rm2 with mi and m2 integers, we can rewrite Eq. 



in the form 



k oxi^v'V I 

Adding to this expression the trivial Cu contribution we get : 



Sq = {N + iV^)-i/2 E ^kAiQ, k) E exp ^((A; - Q)l + k^{l)), (A5) 



k 

where we defined the coherence factor 



sm(A;/2) 
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APPENDIX B: EVALUATION OF THE FERMIONIC CORRELATOR 



The fermionic correlator is given by the expression 

{expikA{l,t) exp-ikA{l',t'))^^, (Bl) 



where the bosonic representation of A{1, t) has been given in Eq. |T6| above. Eq. ^T] may be 
easily evaluated by repeated application of the well-known formula 

exp{A + B) = exp Aexp B exp{-l/2[A, B]), (B2) 

valid for any pair of operators A and B that satisfy [A, [A, B]] = [B, [A, B]] = 0. 



We first use Eq. B2 in Eq. Bl in order to combine the two exponentials into one, 



{expikA{l,t) exp —ikA{l' ,t')) = 
= (exp tk {A{1, t) - A{1\ t')})^^ exp [A{1, t),A{l\ t')]] , 

where the commutator on the right hand side of Eq. |B^ is a c— number given by (cf. Eq. 



T6D : 

[Ail, t), A{1\ t')] = ~2i Y: v\q) sin [t - t')) Re [0(g, /')] • (B4) 



(B5) 



Next, we write A{l,t) — A{1' ,t') as the difference of two bosonic operators, 

A{1, t) - A{1\ t') = B{1, 1'; t, t') + B\l, t, t'), 
B{1, V- 1, t') = v{q) /)e--^* - /')e-'"^*') 

Using once more Eq. we obtain : 

(expiA;(i3 + B^)) = exp -y [i?, B% (B6) 



where 



[B, St] = 2 5: v\q) I - cos M<l\t) Re l)<P{-q, m)] | , (B7) 

and we have used the fact that no bosons are excited in the ground-state. Combination of 
these results leads to the expression quoted in Eq. [18| of the main text. 
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APPENDIX C: TEMPERATURE DEPENDENCE OF THE FERMIONIC 

CORRELATOR 



We show below that the correlation length of the fermion subsystem at finite temperature 
is ^ = vp/T. For this we use a modification of an approach originally due to Lieb, Schultz 
and Mattis It is based on a representation of the instantaneous correlation function 
Cq(/, 0; if: = t') of Eq. |TD| as the Toeplitz determinant 



Cq{1) = {exptQj2H(lJ)pq) = det{M) 



(CI) 



Mi 



kl 



6u - (1 - e''^)gki. 



The Ixl matrix M is constructed out the fermionic correlation functions Qkm = sin(7rc(A; — 
m))/7r(k — m) = {l/2n) jZ-j^dqUqC^''^'^^'"^^ where the Fermi function = 1/(1 + 6^''/"^). Notice 
that here we do not subtract the g = component of the density in Eq. [Cll . Next, we use 
the asymptotic behavior of the above determinant at large Z, 



det{M) -expf/D-fsTs], 



(C2) 



where 21 



D = iQ 



dqdk ^ 



dq ^ 



-n 
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1 — cos kl 



(C3) 



(C4) 



2(27r) 

and iQfiq = ln[l — (1 — e''^)ng]. 

At T = fig = riq and we recover our results in the main text. For finite temperature, 



fir 



is not so simple. For our purposes it is sufficient to expand fig up to terms of order 



g2 

iQUq = iQUq - ^^^9(1 - T^q)i 



(C5) 



from which the corrections to the two terms that appear in the exponent in Eq. |C2| follow : 



ID = iQcl 



(C6) 



dk 1 - cos kl 



Jo k exp{kvF/T) — 1 

sinh7r//e 
= — t: In ; • 

Combining ( |C6D and ( |C7|) , one finds that at Q ~ vr and large / 

|C(0|~exp [-tt/^]. 
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FIGURES 

FIG. 1. The function for the CuO chains for several values of the hole concentration c. 

FIG. 2. Frequency dependence of $(0, CI) for three values of the hole concentration c. The ratio 
of the effective exchange constant to the hole hopping integral is Je///r = 1/3. 

FIG. 3. Equal-time spin-spin correlation of the CuO chains as function of q and of the hole 
concentration c. 
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